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1. Introduction 

The admittance determines the current response of an electrical conduc- 
tor to oscillating voltages applied to its contacts. The admittance and the 
nonlinear dc-transport have in common that they require both an under- 
standing of the charge distribution which is established as the conductor 
is driven away from its equilibrium state. The dynamic and the nonlinear 
transport can thus be viewed as probes of the charge response of the con- 
ductor. In contrast, linear dc-transport is determined by the equilibrium 
charge distribution alone. Thus the dynamic conductance and the nonlin- 
ear transport provide information which cannot be gained from a linear 
dc-measurement . 

It is the purpose of this article to provide a discussion of ac-transport 
and nonlinear transport in small mesoscopic conductors. Our discussion 
leans heavily on the scattering approach to electrical dc-conductance and 
can be viewed as an extension of this approach to treat transport beyond the 
stationary ohmic regime. Our emphasis is on the low-frequency admittance 
for which expressions can be derived which are very general and can be 
applied to a large class of problems. Similarly our discussion emphasizes the 
departure from linear ohmic transport in the case when nonlinearity first 
becomes apparent. The principles which govern the low- frequency behavior 
and the weakly nonlinear behavior apply of course as well for high-frequency 
transport and for strongly nonlinear transport. We underline this fact by 
treating the simple example of a resonant tunneling barrier for its entire 
frequency range and by calculating the nonlinear I-V characteristic for a 
large range of applied voltages. 
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The general guiding principles of our discussion are charge and cur- 
rent conservation and the need to obtain gauge invariant expressions. The 
conservation of charge becomes a fundamental requirement under the fol- 
lowing condition [1]: Suppose that it is possible to find a volume which 
encloses the mesoscopic part of the conductor including a portion of its 
contacts such that all electric field lines remain within this volume. Clearly, 
if the mesoscopic conductor is formed with the help of gates the volume 
must be large enough to enclose also a portion of the gates. The electric 
field is then localized within Such localized electric (and magnetic) field 
distributions are always implicitly assumed when a circuit is represented 
simply in terms of R, C and L elements. If the field is localized in 0, the 
electric flux through the surface of that volume is zero and according to 
the law of Gauss the total charge inside this volume is zero. Consequently, 
application of a time-dependent or stationary voltage to the conductor (or 
a nearby gate) leads only to a redistribution of charge within O but leaves 
the overall charge invariant. Therefore, a basic ingredient of a reasonable 
discussion of transport of an electrical conductor is a theory of the charge 
distribTition with the overall constraint that the total charge is conserved. 
If the total charge is conserved then the currents measured at the con- 
tacts of the sample are also conserved. As we will show, this leads to sum 
rules for the admittance coefficients of a mesoscopic conductor and leads 
to sum rules for the coefficients which govern the nonlinear dc-transport. 
The conservation of the total charge is also connected with the fact that 
an electrical conductor does not change its properties if we shift all poten- 
tials by an equal amount. The principle of gauge invariance states that the 
results can depend only on voltage differences. This leads to a second set 
of sum rules for the admittance coefficients and for the nonlinear transport 
coefficients [2]. 

The three guiding principles which we have taken for the development of 
a theory are thus very closely related and hinge on the localization of electric 
and magnetic fields. The existence of a Gauss volume is not, however, a 
trivial requirement. A different point of view is taken in a recent article 
by Jackson [3]. In his circuit the potential distribution depends on the 
location of the battery vis-a-vis the resistor. Such a point of view would 
make it necessary to develop a theory for the entire circuit including its 
components used to drive it out of equilibrium. We believe that the notion 
of localized fields is a more fruitful one and is closer to the experimental 
reality as it is encountered in mesoscopic physics. In what follows, it is not 
only assumed that the electric fields are localized, but that the localization 
volume fl is in fact of the same dimension as a phase breaking length. This 
permits us to formulate a theory for the mesoscopic structure alone and to 
treat the completion of the circuit as a macroscopic problem. 
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A charge and current conserving approach to mesoscopic ac-conductance 
was developed by one of the authors in collaboration with Pretre and 
Thomas [4]. Similar work, without an effort to achieve a self-consistent 
description, was also presented by Pastawski [5] and Pu and Dudley [6]. 
In Ref. [4] charge and current conservation was achieved by attributing a 
single self-consistent potential to a conductor. A discussion which treats 
the potential landscape in a Hartrcc-like approach was given in Rcf. [7] for 
metallic screening and in Ref. [1] for long range screening appropriate for 
semiconductors. These two discussions can be extended to include an ef- 
fective potential as it occurs in density functional theory [8]. Both the limit 
of a single potential and the case of a continuous potential are important 
but for many applications the first is to crude and the later is to complex. 
In this work we will treat an approach which is in between these two: In 
this generalized discrete potential model the conductor is divided into as 
many regions as is necessary to capture the main features of the charge and 
potential distribution. This permits the description of charge distributions 
which are mainly dipolar or of a more complex higher order multi-polar 
form. The authors applied such an approach to the low-frequency behavior 
of quantized Hall samples [9], quantum point contacts [10] and to nonlinear 
transport [2]. 

We will only briefly review the scattering approach to electrical con- 
ductance, mainly to introduce the notation used in this article. For a more 
extended discussion of some of the basic elements of the scattering approach 
we refer the reader to the introductory chapter of this book or to one of 
the reviews [11, 12, 13, 14]. The advantage of the scattering approach, as 
formulated by Landauer [15], Imry[ll] and one of the authors [16, 17] lies 
in the conceptually simple prescription of how to model an open system, 
i.e. how to couple the sample to external contacts which act as reservoirs of 
charge carriers and provide the source of irreversibility. The self-consistent 
potential in the context of transport was already of interest to Landauer [18] 
and is in fact the main issue which distinguishes his early work from pre- 
ceding works on transmission through tunnel contacts [19]. A more recent 
discussion of the potential distribution in ballistic mesoscopic conductors 
is provided by Levinson [20] and technically is closely related to our ap- 
proach [1]. Por the discrete potential model discussed here, the potential 
distribution can be discussed in a purely algebraic formulation which we 
will present in some detail. A large portion of this article is devoted to an 
explicite illustration of our approach. We discuss the low frequency admit- 
tance of wires with barriers and the low frequency admittance of quantum 
point contacts. In particular, we generalize published work on the low fre- 
quency admittance of the quantum point contact to investigate the effect 
of the gates used to form the quantum point contact. The last section of 
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this work presents a discussion of ac-transport and nonlinear transport over 
a large frequency and voltage regime under the assumption that the only 
important displacement current is that to a nearby gate. 

The discussion of potentials also sheds light on the limits of validity of 
the scattering approach as it is used to describe dc-transport: The scatter- 
ing approach is often classified as a noninteracting theory: That is incorrect. 
Since scattering states can be calculated in an effective potential the range 
of validity of the scattering approach is at least the same as density func- 
tional theory. That makes it possible to include exchange and correlation 
effects and makes the scattering approach a mean-field theory. As is known, 
for instance from the BCS-theory of superconductivity, mean-field theories 
can be very powerful tools which render any remaining deviations very hard 
to detect. 

The low-frequency behavior discussed here should experimentally be 
accessible with much of the same techniques as those that are used for 
dc- measurements. We mention here the work of Chen et al. [21] and Som- 
merfeld ct al. [22] on the magnetic field symmetry of capacitances. We 
mention further a recent experiment by Field et al. [23] who used a quan- 
tum dot capacitivcly coupled to a two-dimensional electron gas to measure 
the density of states at the metal-insulator transition. In contrast, exper- 
iments at high-frequencies require special efforts to couple the signal to 
the mesoscopic conductor and to measure the response. As examples, we 
cite here the work by Pieper and Price on the admittance of an Aharonov- 
Bohm ring [24] , investigations of photo-assisted tunneling by Kouwenhoven 
et al. [25], noise measurements at large frequencies [26], and transport in 
superlattices [27]. We emphasize that already the low- frequency linear ad- 
mittance presents a very interesting characterization of the sample. Indeed, 
we hope, that this article demonstrates that there is considerable room for 
additional theoretical and experimental work in this area. 
Nonlinear dc-effects of interest are asymmetric current- voltage characteris- 
tics and rectification [28] , the evolution of half- integer conductance plateaus 
[29, 30], the breakdown of conductance quantization [31], and negative 
differential conductance and hysteresis [32]. Like in frequency-dependent 
transport most of the theoretical work shows no concern for self-consistent 
effects. An exception is a numerical treatment of a tunneling barrier by 
Kluksdahl ct al. [32]. 

The emphasis on the role of the long-range Coulomb interaction in this 
work should be contrasted with the recent discussions of mesoscopic trans- 
port in the framework of the Luttinger approach. The Luttinger approach 
treats the short-range interactions only. However, since a one-dimensional 
wire cannot screen charges completely, a more realistic treatment should 
include the Coulomb effects. Without the long-range Coulomb interaction 
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the results based on Luttinger models only are not charge and current 
conserving and are not gauge invariant. In this context we note that the 
theory of the Coulomb blockade has been very successful in describing the 
important experimental aspects. The Coulomb blockade is a consequence 
of long-range and not of short-range interactions. In any case the Hartree 
discussion given here is Tiseful as a comparison with other theories which 
include interactions. To distinguish non-Fermi liquid behavior from those 
of a Fermi liquid it will be essential to compare experimental data with a 
reasonable, i. e. self-consistent, Fermi liquid theory. 

2. Theory 

In this section wc recall briefly the theory of ac-transport and weakly non- 
linear transport which has been developed in Refs. [1, 2, 4, 7, 10]. A review 
of these works is provided by Ref. [33]. 



2.1. MESOSCOPIC CONDUCTORS 

Consider a conducting region which is connected via leads and contacts 
a = 1,...,N to N reservoirs of charge carriers as shown in Fig. la. In 
order to include into the formalism the presence of nearby gates (capac- 
itors), some parts of the sample are allowed to be macroscopically large 
and disconnected from other parts. A magnetic field B may be present. It 
is assumed that the leads of the conducting part are so small that carrier 
motion occurs via one-dimensional subbands (channels) m = 1, ...,Ma for 
each contact. Moreover, the distance between these contacts is assumed to 
be short compared to the inelastic scattering length and the phase-breaking 
length, such that transmission of charge carriers from one contact to an- 
other one can be considered to be elastic and phase coherent. A reservoir 
is at thermodynamic equilibrium and can thus be characterized by its tem- 
perature Tq and by the electrochemical potential jia of the carriers. While 
we assume that all reservoirs are at the same temperature T, we consider 
differences of the electrochemical potentials which is usually achieved by 
applying voltages Va to the reservoirs. We fix the voltage scales such that 
Va = corresponds to the equilibrium reference state where all electro- 
chemical potentials of the reservoirs are equal to each other, = (Iq. The 
problem to be solved is now to find the current /^(t) entering the sample at 
contact a in response to a (generally time dependent) voltage Vp{t) at con- 
tact /?. The general time-dependent nonlinear transport problem is highly 
nontrivial and we consider mainly weakly nonlinear dc-transport and linear 
Zow-frequency transport. 



6 



M. BUTTIKER AND T. CHRISTEN 




Figure 1. a) Multi-terminal sample including a gate; the voltages Va induce charges 
{q{x)) which influence the potential {U{x)), and drive currents {Ip). b) Decomposition of 
the local density of states in a two-terminal sample into local partial densities of states. 

2.2. THE SCATTERING APPROACH 

The sample is described by a unitary scattering-matrix S [11, 15, 16, 34]. 
Unitarity together with micro-reversibility imphes that under a reversal 
of the magnetic field the scattering matrix has the symmetry S'^{B) = 
S{—B) [16, 17]. Furthermore, the scattering matrix S for the conductor can 
be arranged such that it is composed of sub-matrices Sap{E, \U{x)\) with 
elements Sa^nmiE, [U{x)]) which relate the out-going current amplitude in 
channel n at contact a to the incident current amplitude in channel m at 
contact /?. The scattering matrix is a function of the energy E of the carrier 
and is a functional of the electric potential U{x) in the conductor [1]. For 
the dc-transport properties, all the physical information which is needed, 
is contained in the scattering matrix. 

The electric potential U{x,{Vy}) is a function of the voltages Vy = {fj,^ — 
jJLQj/e applied at the contacts and at the nearby gates. Thus the scattering 
matrix depends not only on the energy of the scattered carriers but also on 
the voltages, Saf^{E, {Ky}). With the help of the scattering matrix, we can 
find the current dIa{E) in contact a due to carriers incident at contact (5 
in the energy interval (i?, E + dE). It is convenient to introduce a spectral 
conductance [2] 

g^p{E) = ^ TriUE, {V^})6^p - sUe, {V^})s^fs{E, {V,})] (1) 
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such that the current at energy E becomes dIa{E) = ga(}{E){dE/e). The 
unity matrix Iq, Uves in the space of the quantum channels in lead a 
with thresholds below the electrochemical potential. Note that this ma- 
trix is also a (discontinuous) function of energy and of the potential. It 
changes its dimension by one whenever the band bottom of a new sub- 
band passes the electrochemical potential. The current through contact 
a is the sum of all spectral currents weighted by the Fermi functions 
f{z) = (1 -f- exp(2;/A;BT))~^ of the reservoirs [35] at temperature T 

TV 

Ia = J2 (dE/e) f{E - MO - eVp) g^piE, {V^}) . (2) 
13=1'' 

Linear transport is determined by an expansion of Eq. (2) away from the 
equilibrium reference state to linear order in Vy. The linear conductance is 
found to be 

g2] = J dE i-dEfo) gat3{E) , (3) 

where the Gap are evaluated at = ... = V/v = 0. Consequently, the linear 
conductances depend on the equilibrium electrostatic potential U^^'^\x) = 
U{x, {0}) only. In contrast, both the ac-transport and the nonlinear trans- 
port depend explicitly on the nonequilibrium potential. A discussion of the 
nonequilibrium potential requires knowledge of the local charge distribu- 
tion in the conductor. Within linear response, the charge response is related 
to the DOS of the conductor at the Fermi energy. The local DOS can be 
obtained from the scattering matrix by a functional derivative with respect 
to the local potential [1] 



dE Am ^ 

ap 



^'^l^eSUix) e5U{xf''^ 



ap 



As indicated in Eq. (4), the local DOS can be understood as a sum of local 
partial densities of states (PDOS) dn{a,x,P)/dE [36]. The sum is over all 
injecting contacts (3 and all emitting contacts a. The meaning of a local 
PDOS dn{a,x,P)/dE is then obvious: it is the local density of states asso- 
ciated with carriers which are scattered from contact /3 to contact a. The 
local PDOS arc illustrated in Fig. lb. The global DOS, dN/dE, follows by 
a spatial integration of Eq. (4) over the sample region which encloses the 
mesoscopic structure with a boundary deep inside the reservoirs. Integra- 
tion of the local PDOS over the whole sample leads to the global PDOS, 
dNap/dE. Clearly, it holds dN/dE = Y.a(^ dN^fs/dE. Note that Eq. (4) 
counts only scattering states. Pure bound states, e.g. trapped or pinned at 
an impurity, are not included. 



8 



M. BUTTIKER AND T. CHRISTEN 



The local PDOS dn{x,P)/dE which contains information only on the con- 
tact from which the carriers enter the conductor (irrespectively of the con- 
tact through which the carriers exit) is called the injectivity of contact 
(3 and is given by dn{x,j3)/dE = '^^dn{a,x,l3)/dE. The local PDOS 
dn{a, x)/dE which contains information only on the contact through which 
carriers are leaving the sample but contains no information on the contact 
through which carriers entered the conductor is called the emissivity of 
contact a and is given by dn{a,x)/dE = Y^pdn{a,x, (5) / dE. Due to reci- 
procity, the injectivity at a magnetic field B is equal to the emissivity at 
the reversed magnetic field, dnsix, a)/dE = dn-B{ci,x)/dE. 
Throughout this work we use the discrete potential approximation in which 
the conductor is divided into M regions ^Ik, k = 1,...,M. The charge 
and the potential in region k are denoted by and U^- Since we are 
are interested in the charge variation in response to a variation in volt- 
age, we introduce the DOS of O^. For later convenience, we express these 
DOS in units of a capacitance by multiplication with e^. Thus, = 
j^^dx {dn{x)/dE) is the DOS of region k. The PDOS of region k are 
E)ak(3 = /qj. dx{dn{a, x, (3) /dE). Obviously, the injectivities and the emis- 
sivities are D^p = /f2^ dx{dn{x, 0)/dE) and Dak = In^. dx{dn{a,x)/dE), 
respectively. 

2.3. LINEAR LOW-FREQUENCY TRANSPORT 

We are interested in the admittance Gajsioj) which determines the Fourier 
amplitudes of the current, 6Iaii^), at a contact a in response to an oscil- 
lating voltage 6V/3{u;) ex.p(—iLot) at contact /3 



To investigate the low-frequency limit, we expand the admittance in powers 
of frequency 



The dc-conductance G^^ has already been derived in Eq. (3) . The first order 
term is determined by the emittance matrix Ea^ = i{dGap/duj)^=Q. The 
emittance governs the displacement currents of the mesoscopic structure. 
The second order term, Ka/^, contains information on the charge relaxation 
of the system. As an example, consider for a moment a macroscopic capac- 
itor C in series with a resistor R. For this purely capacitive structure with 
vanishing dc-conductance, G = 0, the emittance is the capacitance, E = C, 




(5) 



(6) 
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and K = R^'^^C'^ contains the RC time with the charge relaxation resis- 
tance i?^^) = R. More generally, for any structure consisting of N metallic 
conductors, each connected to a single reservoir, the emittancc matrix is 
just the capacitance matrix. These simple results must be modified for 
mesoscopic conductors and conductors which connect different reservoirs. 
Firstly, it is not the geometrical capacitance but rather the DOS dependent 
electro-chemical capacitance which relates charges on mesoscopic conduc- 
tors to voltage variations in the reservoirs [37]. Secondly, conductors which 
connect different reservoirs allow a transmission of charge which leads to 
inductance-like contributions to the emittance [1, 4, 10]. Thirdly, the charge 
relaxation resistance cannot, in general, be calculated like a dc-resistance 
[37, 38]. 

To illustrate our approach, we derive now a general expression for the 

emittancc. Wc first notice that for the purely capacitivc case the cur- 
rent (displacement current) at contact a is the time derivative of the to- 
tal charge SQa on the capacitor plates connected to this contact. Hence, 
5Ia = —iuj5Qa and the emittance is given by = dQa/dVp. If trans- 
mission between different reservoirs is allowed, the charge 5Qa scattered 
through a contact can no longer be associated with a unique spatial region. 
The charge at a given point is rather injected at different contacts and 
ejected at different contacts. However, the charge emitted at a contact can 
still be calculated within the scattering approach. We decompose it in a 
bare and a screening part 6Qa = SQa^ + 6Qa ^ . The bare part of the charge 
corresponds to the charge which is scattered through the contact a for fixed 
electric potential and is thus given by 

dQ(^^ = ^D^p6V^ , (7) 

where Dafs = J2k ^akfi is the global PDOS. The screening part, on the 
other hand, is associated with the charge which is scattered through con- 
tact a for a variation in the electric potential only. Since a shift of the 
band bottom contributes with a negative sign and since the potential is 
in general spatially varying, SQa'' is connected to the local PDOS by 

is) 

5Qa = —J2k-^ctk^^kj where D^k is the emissivity of region k into con- 
tact a. If we introduce the vector of emissivities = {Dai, ■■■,DaM) and 
write the potential variation as a vector, 5\J = {6Ui, ....,6Um), the charge 
emitted into contact a is 



(8) 
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Figure 2. Change (dashed) of the equihbrium potential (sohd) in a two-terminal sample 
due to a voltage variation 5V\ (dotted) in the left reservoir. 



In linear response the potential variation 6\J is linearly related to the po- 
tential variations 5Vp at the contacts of the sample. We write 



<5U = ^u, 5V^ 



(9) 



where the response coefficients Uq, = {uia, ■■■■,UMa) are called the charac- 
teristic potentials of contact a. The emittance can be written as 



Eal3 — Daf3 - DqU/j 



(10) 



To complete the calculation of the emittance, we next need a discussion of 
the characteristic potentials. 



2.4. PILED-UP CHARGE AND SELF-CONSISTENT POTENTIAL 

The variation of the charge iJg^ is the sum of a bare charge and a screening 
charge. The bare charge can be expressed with the help of the injectivity 



(11) 



The screening charge is connected to the electrostatic potential [20] via the 
Lindhard polarization function, which is here a matrix 11 with elements 



(^q(«) = -n (5U 



(12) 
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The Lindhard matrix can be expressed in terms of the scattering states. 
It is in general not a diagonal matrix, i.e. the charge response is in gen- 
eral nonlocal. Non-local effects are, however, small quantum mechanical 
effects which can be neglected deep inside a conductor with a sufficiently 
large electron density. In a quasiclassical treatment the Lindhard matrix 
becomes local, 11^; = d^iDk- Here is the local DOS in region k. 
The total charge, (5q = Jq*-*^ + ^q*-''-' , acts now as the source of the nonequi- 
librium electric potential in the Poisson equation. For our discrete model 
we also have to discretize the Poisson equation. We introduce a geomet- 
rical capacitance matrix C which relates the charges to the electrostatic 
potentials, (5q = C6\J. If one uses this matrix C and expresses the electric 
potential in terms of the characteristic potentials, one can write the Poisson 
equation in the form (C-l-n)ua = D^^ which determines the characteristic 
potentials Uq, 

u, = (n + C)-iD« . (13) 
Going back to Eq. (10) we find for the emittance 

E^f, = D^^-B^^\n + C)-'nf . (14) 

We notice that the bare contribution to the emittance occurs with a positive 
sign and the Coulomb induced contribution occurs with a negative sign. 
Depending on which contribution dominates we call the emittance element 
capacitive {E^a > 0, Eai3 < for a 7^ /3) or inductive-like {E^a < 0, 
Eap > for a / /3). 

The nonequilibrium charge distribution becomes 

(5q = C5U = ^D^'^JF^ -roU . (15) 



This determines the nonequilibrium steady-state to linear order in the ap- 
plied voltages. The nonequilibrium charge-distribution can be written in 
terms of an electrochemical capacitance matrix C^,fe/3 [2, 9, 10] which de- 
termines the net charge variation in region k in response to a potential 
variation at contact /3. In vector notation, we have 

6ci = ^C^,^SVfs , (16) 

where 

C^,;3 = Cu;3 = C(n + C)-^D^') . (17) 

We do not present here a general discussion of the charge relaxtion term 
Kap. This requires a dynamical screening theory, and we have results only 
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for special cases. We discuss one of these special cases in Sect. 2.7 and Sect. 
3.3.2. 

2.5. WEAKLY NONLINEAR DC-TRANSPORT 

The nonequilibrium potential distribution is not only of importance in ac- 
transport but also in nonlinear transport. Knowledge of the nonequilibrium 
potential distribution to first order in the applied voltages permits to find 
the nonlinear I-V-characteristic up to quadratic order in the voltages, 

= E G%Vp + G^X^pV, + OiV') . (18) 

The coefficients G^'^^ and G^^J]^ are obtained from an expansion of Eq. (2) 
with respect to the voltages Va- One obtains for the linear conductance Eq. 
(3). An expansion of Eq. (2) up to 0{V'^) yields G^ = (1/2) / dE {-dEf) 

{'^dv^9af3 + edEgaisSis^) ■ Writing dv^ga/s in terms of the derivatives dga/s/dUk 
and the characteristic potentials yields 

= \jdE {-dEf) E ^ (2n., - 5,,) . (19) 

Expressing the characteristic potentials in terms of the injectivities we find 
that the rectification coefficient is given by 

= \j dE{-dEf) (2(Vui?a/3)* (n + C)-iD» + eSp^dEgap) . (20) 

For a quantum point contact, Eq. (20) has been discussed in Ref. [2]. In 
Sect. 3.3 we will apply this result to the resonant tunneling barrier. 

2.6. CHARGE CONSERVATION, GAUGE INVARIANCE, AND MAGNETIC 
FIELD SYMMETRY 

Since the system under consideration includes all conductors and nearby 
gates, the theory must satisfy charge (and current) conservation and gauge 
invariance. Due to micro-reversibility the linear response matrices must 
additionally satisfy the Onsager-Casimir symmetry relations. 
Let us first discuss charge conservation, which states that the total charge 
in the sample remains constant under a bias. This implies also current 
conservation, J2a = 0- Imagine a volume J7 which encloses the entire 
conductor including a portion of the reservoirs which is so large that at 
the place were the surface of J7 intersects the reservoir all the characteristic 
potentials are either zero or unity. According to the law of Gauss, one 
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concludes that the total charge remains constant. Application of a bias 
voltage results only in a redistribution of the charge. If the conductor is 
poor, i.e. nearly an insulator, the contacts act like plates of capacitors. In 
this case long-range fields exist which run from one reservoir to the other 
and from a reservoir to a portion of the conductor. But if we chose the 
volume Q to be large enough then all field lines stay within this volume. 
Charge and current conservation imply for the response coefficients the sum 
rules 

= E^S = = E^a/3 = Y.G%. = • (21) 

k cx cx oc a 

Second consider the fact that only voltage differences are physically mean- 
ingful. Gauge invariance means that measurable quantities are invariant 
under a global voltage shift SV in the reservoirs, 5Va ^ ^Va + SV. A global 
voltage shift corresponds, of course, only to a change of the global voltage 
scale. Consequently, the characteristic potentials satisfy [37] 

Eu° = l ■ (22) 

a 

For the response coefficients the following sum rules hold 
E C,,kp = E = E = E = E(gS, + = . (23) 

(3 f} 13 

Note that due to charge conservation and gauge invariance, the geometrical 
capacitance matrix C has the zero mode 1 = (1, 1, ....1), corresponding to 
a unit potential in each region. Thus C cannot be inverted. However, the 
Green's function (11 -|- C)~^ which solves the Poisson equation exists (see 
Eq. (13)). Expressing the characteristic potentials with the help of the 
Green's function and the injectivities gives 

l = (n + C)-^D (24) 

where D = (Di, D2, Dm) is the vector of the local DOS. Equation (24) 
is the key property of the Green's function needed to show that our final 
results are charge and current conserving. From Eq. (24) it follows imme- 
diately that the sum over the elements of a row of the Lindhard matrix is 
equal to the DOS in region k, i.e. Dk = J2i ^ki- Since the Lindhard matrix 
is symmetric, it holds also Dj. = ^i^ik- The above mentioned statement, 
that in the quasiclassical local case 11^/ = Dkdki holds, is now clear. 
In the presence of a magnetic field, the admittance matrix must additionally 
satisfy the reciprocity relations 



(25) 
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which is a consequence of time-reversal symmetry of the microscopic equa- 
tions and the fact that we consider transmission from one reservoir to an- 
other [1, 16, 17]. These Onsager-Casimir relations are only valid for the 
linear response coefficients, i.e. close to equilibrium. Of course, the symme- 
try relations (25) hold individually for G^'J], Eap, and Kap- Note that the 
emittance is symmetric, Eap = Epa, in the purely capacitive case. 

2.7. FREQUENCY-DEPENDENT TRANSPORT: SINGLE POTENTIAL 
APPROXIMATION 

In this section we show that it is possible to discuss the full frequency- 
dependence of the admittance, if the electric potential in the mesoscopic 
conductor can be approximated by a single variable. The external response 
is defined as the response for fixed electrostatic potential (i.e. the 'bare' 
response). A general result for the external response has been derived in 
Ref. [4]. If the differences of wave vectors k{E) and k{E + hu) deep in the 
reservoirs are neglected, the external response can be expressed in terms of 
the scattering matrix only, 



The scattering matrices and the Fermi functions are here evaluated at the 
equilibrium reference state. The real part of the ac-conductance, (26) is 
related via the fluctuation-dissipation theorem to the current-current fluc- 
tuation spectra derived in Ref. [39]. An expansion of Eq. (26) to linear order 
in CO gives [4] 



where the dNaf^/dE are the global PDOS, in which the functional derivative 
with respect to the local potential (and the integration over the volume O) 
is replaced by a (negative) energy derivative. Asymptotically, for a large vol- 
ume the global PDOS obtained from an energy derivative is identical to 
the integral of the local PDOS associated with a functional derivative with 
respect to the potential according to Eq. (4). However, for a finite volume, 
there are typically small differences of the order of the Fermi wavelength 
divided by the size of the volume $1. For a more detailed discussion, we 
refer the reader to Ref. [36]. 

In order to obtain the full admittance, we have still to find the internal 
response (the 'screening' part of the admittance). A general result for the 
internal response is not known. In the simple case, where the conductor 
is treated with a single discrete region (M = 1), the external response 



Glp{o^) = dETr[laSap - slgiE)s^f^{E + hcj)] 



{f{E)-f{E + huj)) 



(26) 




(27) 
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defines, however, also the internal response. We assume that the sample 
is in close proximity to a gate which couples capacitively to the conduc- 
tor with a capacitance coefficient C. The current response at contact a is 
5Ia{oo) = X] G^^(u;)(5Va + GQ(a;)(5C/i where G]^{(jo) is the (unknown) internal 
response of the conductor generated by the oscillating electrostatic poten- 
tial 5U\. The current induced into the gate is 5Ig{(jo) = —iujC{5Vg — 5Ui). 
Now we can determine G'q(<^) from the requirement of gauge invariance. In 
particular, if all potentials are shifted by —5Ui, it follows immediately that 
G^ail^) = ~ ^a/jC*^)- Using current conservation, X^^, 5Ia + 5Ig = 0, we 
find for the conductances [4] of the interacting system (a,/3,7, (5 ^ g), 

^^/^^"^ - -l + i^/.C)E,sG^^,i.y 

""-^"^ - "i + (V-C)E,,G^,(-)' ^''^ 
GUu^) = ...fi^fS^^l .. . (31) 



l + {i/ujC)j:-ysG'^si^ 



The single-potential approximation provides a reasonable description only if 

the conductor has a well-defined interior region, which might be described 
by a single uniform potential. Furthermore, it is assumed that the only 
relevant capacitance is that of the sample and the gate. Examples for which 
this approach is reasonable are quantum wells or quantum dots or cavities 
for which the Coulomb interaction is so weak that single electron effects can 
be neglected. In Section 3.3 we discuss the case of a double barrier with a 
well which is capacitively coupled to a gate. 



3. Examples 

In a first example we discuss the omittance of a wire which contains a single 
barrier. The wire is capacitively coupled to a gate. As a second example we 
discuss the capacitance and emittance of a quantum point contact formed 
with the help of gates [10]. It turns out that steps in the capacitance and 
the emittance occur in synchronism with the well-known conductance steps. 
As a third example, we treat the resonant tunneling barrier close to a res- 
onance, for which we calculate the nonlinear current-voltage characteristic 
and the frequency-dependent conductance [2] with the single-potential ap- 
proximation. For all examples, we use explicit expressions for the Lindard 
function (matrix) 11 given within the local Thomas-Fermi approximation. 
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gate 




Figure 3. Quantum wire containing an impurity and coupled to a gate. 

where 11 is diagonal with elements determined by the local DOS. In all 
examples we consider finally the zero-temperature limit. 

3.1. ADMITTANCE OF A WIRE WITH A BARRIER 

Our first example is a mesoscopic wire which connects two contacts and is 
capacitively coupled to a nearby gate (see Fig. 3). The wire is assumed to 
be perfect with a uniform potential along its main direction, except for a 
single scatterer, a barrier, or an impurity. The discussion is carried out in 
the semiclassical limit where Friedel oscillations are neglected. We take the 
potential of the impurity to be of short-range and assume that the potential 
to left and right of the scatter can be assumed to be uniform. We assume 
here that the gate couples to the wire only and neglect scattering at the 
interface of the wire and the reservoir. Furthermore, to simplify the discus- 
sion we assume in this section that the capacitance between the charges to 
the left and the right of the barrier can be neglected (i.e. vanishingly small 
C12). Consider first a perfect wire (T = 1) of length L. The gate couples 
with a geometrical capacitance c per length. The DOS per length for a 
one-dimensional wire is 

dn/dE = 2/hv (32) 

where v = (2/m)^/^(£' — eUY/"^ is the velocity of carriers with energy E. If 
the wire is filled up to the electrochemical potential /U of the reservoirs the 
total number n of states per length is 

n{^l) = r dE{dn/dE) = e-\2/TT^aB)^^Hfi - eUf^. (33) 

JeU 
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Here we have introduced the Bohr radius ub = /me^. To determine this 
density self-consistently we consider the Poisson equation. The difference 
of the electron charge density en and the ionic background charge density 
—en"'" is equal to the capacitive charge, 

en{n) - en+ = c{U - Ug) , (34) 

where Ug is the electrostatic potential of the gate. Similarly, near its surface, 
the deviation of the electronic gate charge away from charge neutrality is 
equal to —c{U — Ug). Now we take the gate to be a macroscopic conductor. 
Therefore, the electrostatic potential and the electrochemical potential are 
everywhere related by = Epg + eUg, where Epg is the chemical potential 
(Fermi energy) of the gate. We eliminate the gate potential Ug in Eq. (34) 
and find 

en{ij) = {c/ e){e^n'^ /c + Epg — ^g + eU) . (35) 

Measuring energies from the band bottom of the one-dimensional wire and 
setting eU = gives the electrochemical potential as a function of the gate 
potential, 

^x = {c/ef{T^^aB/2){e^n+/c + EFg-ligf . (36) 

In the limit c ^ the electrochemical potential of the wire is simply de- 
termined by charge neutrality, fx = {tt'^qb (n"*")^. For fJ-g = fJ-o = 
e^ri^ /c + Epg the electrochemical potential vanishes and the wire is empty 
for fjLg > ^Q. The DOS of left moving carriers, drii/dE = l/hv, is di- 
rectly related to the phase accumulated by the electrons traversing the wire, 
Ldni/dE = L/hv = {l/Tr)d(j)/dE. The phase accumulated by a traversing 
electron (at the Fermi energy) is thus given by 

(f) = TrL{c/e^){e^n+/c + EFg- fXg) . (37) 

The scattering matrix element for the transmitting channel is given by 
S21 = cxp{i(t)) which shows that the scattering matrix and thermodynamics 
are intimately related. At the Fermi energy of the wire, the DOS is given 

by 

dn/dE = {l/2'K'^aB){l/c){e^n' /c + Epg- ^lgy^. (38) 

The true electron density in a narrow energy interval at the Fermi energy 
is found by considering a small variation of the electrochemical potentials 
(5/x and 5 jig and the electric potential 5U away from a reference state (//*^^^) , 

^JLg'^\ U^^"^^) which obeys Eq. (35). Variation of this equation gives 



6q = e{dn/dE){6n - e6U) = c{6U - d/ig/e). (39) 
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Solving for 6U gives for the charge 6q = Cg{dn/e — d/ig/e) with an electro- 
chemical wire-to-gate capacitance 

Cg^ = c-^ + {e\dn/dE))-'^ . (40) 

Here the DOS appears in scries with the geometrical capacitance It depends 
on the reference state and is given by Eq. (38) with fj,g = fx'g'. 
Next, we discuss a wire containing a symmetric barrier. At equilibrium, 
and if wc can treat the potential created by the barrier as short-range, the 
DOS to the left and right of the barrier remains unchanged. Hence the 
consideration given above for the electrostatic equilibrium potential also 
characterizes the wire to the left and right of the barrier. Similarly, for a long 
wire the gate-to-wire capacitance is essentially unaffected by the presence 
of the barrier. In the case of transport, on the other hand, the chemical 
potentials of the reservoirs connected to the wire differ by a small amount 
implying a variation of the charge distribution. We introduce two regions Ui 
and of size L/2 to the left and the right of the barrier, respectively. With 
the help of Fig. lb, the semiclassical local PDOS, Dakj}, can be constructed 
with simple arguments. For example, D211 is given by the transmission 
probability times the DOS of VLi associated with carriers with positive 
velocity, hence -D211 = TDi/2.^ with Di = D2 = D/2 due to symmetry. 
To determine D212, we note that in the semiclassical limit considered here, 
there are no states in Jli associated with scattering from contact 2 back to 
contact 2, hence it holds D212 = 0. With similar arguments one finds for 
the semi-classical PDOS 

Dakp = Dk (r/2 + 5^p{R 5^k - T/2) ) , if a, /3 7^ 3 . (41) 

Form Eq. (41) we obtain for the emissivities and injectivities Dli = D\i = 
+R)D and Df^ = = (1/4)TD. The integrated geometrical ca- 
pacitance is denoted by C = cL. With these quantities the Poisson equation 
reads 

{D/2){l + R)6Vi + {TD/2)5V2-D5Ui = C{5Ui - SUs) 

{T D /2)5Vi + {I + R){D / 2)5V2 - D5U2 = C{5U2-5U^) (42) 

These equations can also be derived by noticing that the variations SF^ 
of the quasi-Fermi levels [34] in the regions k = 1 and 2 are given by 
6F1 = {{l+R)6Vi+T6V2)/2 and5F2 = {il+R)dV2+T6Vi)/2, respectively. 
For a local Lindhard function, the charge in region k is then Sqk = D^.{SFk — 
5Uk) = {C/2){5Uk-5U^). 

For the charge neutral case, C = 0, Eqs. (42) are familiar from Landauer's 
discussion of the potential drop across a barrier [15, 18, 34]. We assume 



ADMITTANCE AND NONLINEAR TRANSPORT 



19 



again that the gate is a massive conductor for which we have 113 = Epg+ells 
everywhere. We use Eqs. (42) to evaluate the characteristic potentials and 
get 

li_i + R)D 1 TP C 

Due to symmetry, it holds U22 = U2\ = u\2 and U2z = uia- Defining 
Ay = 5V\ — 5V2, the potential drop across the barrier can be written in 
the form 



MJ = SUi - 6U2 = {Cg/C)RAV (44) 

where Cg = Lcg is the electro-chemical gate-to-wire capacitance given by 
Eq. (40). The voltage difference AU is proportional to the reflection prob- 
ability R and the chemical potential difference of the reservoirs. The total 
charges to the left or to the right of the barrier are the sum of the injected 
charges and the induced charges 

dqk = J2iDki3-DkUkp)6Vp . (45) 

We find, e.g., for 5qi, 

Sqi = Cg{{l/2){1 + R)6Vi + {T/2)6V2 - 6V3}. (46) 

The difference in the charge density on the left and the right hand side 
becomes 



Aq = Sqi - 6q2 = CgRAV. (47) 

In the limit of vanishing gate-to-wire capacitance, the charge difference 
vanishes also. Away from the barrier the wire is charge neutral. In this 
charge-neutral limit the potential difference is determined by the reflection 
coefficient only, 5U = RAV . 

Prom Eq. (46) we find the electrochemical capacitance coefficients C^^ka- 
We write this capacitance matrix in the form 





( c^ 


Cg - C^ 


-Cg 




Cg - C^ 




-Cg 




K -Cg 


-Cg 


2Cg 



where 



C^ = (l/2)(l + i?)Cg. 



(49) 
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TABLE 1. Gate capacitance and capacitance, emittance, polarization, 
and voltage drop of the quantum point contact with a grounded gate for 
limiting cases 









Co ^0 




Co = const. 


C = const. 





Co 





D/2 


l/{{D/2)-'+C-') 


Cm 




(D/4)(l + i?) 


(1 + R)Cg/2 


Ell 


RC^, - DT'/A 


RD/2 


{CJC){RC - DT'/A) 


Aq/AV 


2Cp 


RD/2 


RCg 


AU/AV 


C^/Co 





RCg/C 



It follows C^,i2 = (l/2)TCp. Using the characteristic potentials and the 
PDOS given by Eq. (41) we can calculate the emittance, 

^11 = E22 = {CglC){RC-DT''/A) 
Ei2 = E2i = {Cg/C){TC + DTy4) 
Ei3 = Esi = E23 = E32 = —-£'33/2 = —Cg (50) 

If the transmission probability vanishes, the emittance matrix is purely ca- 
pacitive. In the case where the geometric gate capacitance C tends to zero, 
the electrochemical gate capacitance and the electrochemical capacitance 
across the barrier vanish Cg = = Q but the ratio Cg/C tends to 1. The 
wire is then charge neutral and the charge imbalance Ag vanishes. In the 
single-channel case considered here, the voltage drop is A?7 = RISV . In 
this case the wire responds inductively for all values of the transmission 
probability E = —T^D/A. In the absence of a barrier (R = 0) the wire has 
also inductive emittances E = En = E22 = —E12 = —E21 = —D/A. The 
results of this section are summarized in the rightmost column of table 1. 

3.2. EMITTANCE OF A QUANTUM POINT CONTACT IN THE 
PRESENCE OF GATES 

A quantum point contact (QPC) is a small constriction in a two-dimensional 
electron gas which allows the transmission of only a few conducting chan- 
nels [40, 41]. We consider a symmetric QPC with two gates as shown in Fig. 
4a. and ask for its capacitance and low-frequency admittance. To discretize 
the QPC, we define to the left and to the right of the constriction two 
regions fii and VL2 with sizes of the order of the screening length. If only a 
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Figure 4- a) Quantum point contact with gates. The regions fli^2 at the left and right 
of the point contact can be charged, b) Capacitance Cfi and emittance En of the QPC. 
Capacitances and emittances are in units of femto farads. 



few channels are open, the equihbrium potential of the constriction has the 
form of a saddle [42]. Away from the saddle towards the two-dimensional 
electron gas, the potential has the form of a valley which rapidly deepens 
and widens. In contrast to the previous example, l^i and are now re- 
gions immediatley to the left and right of the saddle, where the equilibrium 
potential and the equilibrium density are not uniform. 
Consider a symmetric QPC. The two gates are taken to be at the same 
voltage VJj. Thus in effect, the two gates act like a single gate. As in the 
previous section, the gate will be treated as a macroscopic conductor. For 
the QPC contact in the presence of the gates the charge distribution is 
not dipolar. It is the sum of the charges in Oi, and at the gates which 
vanishes, 6qi + 6q2 + Sq^ = 0. Assume a single open channel. Due to the 
symmetry of the problem, the geometric capacitance matrix can be written 
in the form 



Co + C 


-Co 


-c 


-Co 


Co + C 


-c 


-c 


-c 


2C 



where Co is the geometric capacitance between ili and 0,2 and where C is 
the geometric capacitance between these regions and the gates. The elec- 
trochemical capacitance has again the form (48), with equal to the 
electrochemical capacitance across the QPC and with Cg being the electro- 
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chemical gate capacitance. With the help of Eq. (17) and the characteristic 
potentials Ukp = (-Dfe/? - Cfj_^kp)/Dk one finds 

= C-i + (D/2)-i ^^^^ 

_ RCo + {C/2){l + E) + 2C,C,D-' 

~ 1 + 2(2C7o + • ^ ^ 

The emittance elements (14) become 

^11 = = - £»t2/4 + CgT/2 

E\2 = = Cg — Ell 

Ei3 = E^i = E23 = E32 = — -E'33/2 = —Cg (54) 

with D = 2Di. The charge difference Aq and the electrostatic voltage drop 
AU = 6Ui — SU2 across the QPC are, respectively, given by 

Aq_ ^ DR{2Co + C) 

AV D + 2i2Co + C) ' ^ ^ 

— = (^6) 
AV D + 2{2Co + C) ' ^ ' 

The discussion in the previous section is a limiting case of the general re- 
sults given here. In the limit Co — > we obtain the results discussed for 
the quantum wire with a barrier collected in column three of table 1. In 
the limit C — > overall charging of the QPC is prohibitive. The charge 
distribution is dipolar and the results of column 1 of table 1 apply. In the 
limit C ^ 00 (corresponding to a geometrical capacitance of the QPC to 
the gates which is much larger than the capacitance across the QPC) an 
electrostatic voltage difference across the QPC does not develop, even so a 
charge imbalance exists. The results for this limiting case are collected in 
the middle column of table 1. 

To generalize these results to the many channel case, we mention that for 
the potential considered the contributions of each individual channel just 
add. Each occupied subband j contributes with a transmission probabil- 
ity T^-'-' to the total transmission function T (= Tn) = '^jT^^\ The dc- 

conductancc (3) is then given by cf} = = -Cf^ = -G^^i = (2eV/i)'^- 
Clearly, G^^^^ vanishes whenever one of the indices corresponds to the gate 
index 3. Similarly the local PDOS (41) d'^^^ of each individual channel sim- 
ply add. Hence, the local DOS becomes = J2j D^k \ and T = 1—R is now 
an average transmission probability defined by T = = D7^ Y,, d9^T^^^ 
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where (k = 1,2). Note that the average transmission probabihty T (7^ T) 
has nothing to do with the dc-conductance. If only a few channels are open 
the potential of the QPC has the from of a saddle point [42]. As a spe- 
cific example we consider a quadratic potential U{x) = Uo{b'^ — x'^)/b'^ if 
|a;| < b, and U{x) = if 6 < |x| < Z. The transmission probabilities T^^^ and 
the PDOS can be calculated analytically from a semiclassical analysis [43]. 
For simplicity, we assume a constant electrostatic capacitance Co = 1 fF 
between Qi and 0,2 and a fixed number of occupied channels in these re- 
gions. Both the height Uq and the spatial variation of the potential should 
in principle be found by minimizing the grandcanonical potential for the 
system [44]. Here we consider Uq as an independent control parameter. We 
assume that no additional channels enter into the regions during the 
variation of Uq. In Fig. 4b. we show the results for a constriction with 
b = 500 nm, I = 550 nm, and with three equidistant channels separated by 
Ep/3 = 7/3 meV. The dotted curve represents the transmission function 
which determines the dc conductance. At each step a channel is closed. 
The dashed and solid curves correspond to a gate with coupling C = 
and C = Cq. respectively. The curves represent the capacitance and 
the cniittancc En. For the two-terminal QPC (C = 0) the capacitance 
vanishes and the emittance is negative for small Uq where all channels are 
open. At each conductance step, the capacitance and the emittance increase 
and eventually merge when all channels are closed. Due to a weak logarith- 
mic divergence of the WKB density of states at particle energies E = eUo 
(where WKB is not appropriate) , the emittance shows steep edges between 
the steps. In the presence of the gate (C 7^ 0), the curves are shifted up- 
wards due to a capacitive contribution of the gate. 



3.3. THE RESONANT TUNNELING BARRIER 



In the framework of the single-potential approximation introduced in Sect. 
2.7 we discuss now the fully nonlinear I-V characteristic and the fully 
frequency-dependent admittance of a resonant tunneling barrier (RTB) 
with a single resonant level. For a comparison of the results of our Hartree- 
like discussion (which to be realistic needs to be extended to a many level, 
many channel RTB) with theories that treat single electron effects, we 
refer the reader to the works by Brudcr and Scholler [45], Hettler and 
Scholler [46], Stafford and Wingreen [47] and Stafford [48]. Consider the 
RTB sketched in Fig. 5a, where two reservoirs are coupled by a double bar- 
rier structure. The well between the barriers is considered to be the single 
relevant region. Additionally, the well may couple electrically with a capac- 
itance C to an external gate. We label quantities associated with the well 
with an index 0. The band bottom of the well between the barriers and the 
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energy of the resonant level are denoted by eUo and by £'0 = eUo + E^, 
respectively. The electrochemical potentials of the particles which are scat- 
tered at the RTB arc assumed to be close to the resonant level, such that 
the energy dependence of the scattering properties is described by a Breit- 
Wigner resonance with a width V. The asymmetry of the barrier is denoted 
by Ar = Fi — where Ta/fi is the escape rate of a particle trapped in 
the well through the barrier a . The scattering matrix elements have then 
a pole associated with a denominator A{E) = E — Eq + ir/2, 



-a/^ = - ^^^J . (57) 

Here, the ipa are arbitrary phases. We assume that the only energy depen- 
dence occurs in the resonant denominator, while the Fq, and the (^q, are 
energy independent. One finds a transmission probability 

TiE) = l-R{E)^^, (58) 

with a maximum value Tmax = 4rir2/r^. The injectivities and emissivities 
are equal to each other due to the absence of a magnetic field and are 



i^0a(^)=i^.0(^) = |^^ . (59) 



Summing up injectivities (or emissivities) yields the total DOS in the well 

Do = e^r/27r|A(£;)p . (60) 
Again, the (P)DOS are here represented in units of a capacitance. 
3.3.1. Nonlinear transport 

Let us discuss the nonlinear I-V characteristic [2] for the charge neutral 
case (C = 0). Note that only the effect of the resonant level is considered 
and that transport is still phase coherent. First, we determine the self- 
consistent potential. The relation between the electrostatic potential shift 
AC/q = Uq — ulf^^ and the voltage shifts Va = (Ma — /^^^'^)/e is obtained 
from the charge-neutrality condition 

/ DoiiE,Uo)dE+ Do2{E,Uo)dE- DoiE,U^''^) dE = , 

J— 00 J —00 J— 00 

(61) 
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Figure 5. a) Resonant tunneling barrier with a single level and biased by voltages 
Vi = {fii — fi^'"'^)/e and V2 = 0. The well may be coupled capacitively with a capacitance 
C to a gate (dark region), b) Nonlinear current-voltage characteristic with Fi = 0.5 meV, 
r2 = 1 meV, AE = (solid), AE = -F (dashed), A^; = -2r (dotted). Thin dotted 
lines indicate the quadratic approximation. 



which can be integrated analytically. The result can be written as a function 
which determines implicitly W = AUq — (Vi + V2)/2 as a function of AV 
only, 

AE, ^ ,AE - e(W - AV/2) 
r arctan[^;y^J = Fi arctan[ ^^y^ 



,AE-e(W + AV/2)^ 
T2 arctan[ — J 



r/2 

Equation (2) yields for the current I = h = —I2 

2eVxT2( ,AE-e^ - AVn), 

' arctan 



(62) 



_ lllM^fl™^)') , (63) 

where AE = /i*-*^^) — {elllf^^ + E^) is the equilibrium distance between 
the Fermi energy and the resonance [2]. Without taking into account the 
self-consistent shift AUq one would get a wrong result which is not gauge 
invariant. The current given by Eq. (63) saturates at a maximum value 
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proportional to 7r/2 — arctan(2A£^/r). The conduction is optimal for AE = 
and Fi = r2 when / = (e//i)rarctan(eAy/r). In Fig. 5b we have plotted 
the characteristic for an asymmetry AF /F = — 1/3 and for various values of 
AE. Due to the complete screening, the resonant level floates up or down 
to keep the charge in the well constant. An expansion of the current yields 
[2] gS?i = -{e^/h){AT/2r)dET (thin dotted lines in Fig. 5b). The case of 
incomplete screening can similarly be treated with our approach. At large 
voltages, the resonance can then eventually fall below the conductance band 
bottom of the injecting reservoir as is known from semiconductor double- 
barrier structures. Finally we mention that, in general, even an clastically 
symmetric resonance can be rectifying if screening is asymmetric. 

3.3.2. AC-Response 

In order to calculate the admittances (28) we must know the external re- 
sponse Gl^p defined for fixed potential. Using the specific scattering matrix 
elements (57), Eq. (26) gives 

G^f^ico) = ijdE m-fiE+n^) (s^^ _ si^iE)s^^{E + hu)) 

= a-ap 91 {(^) (64) 



with 



'^"^^'"^ = rx^^r^ I -^riF2 F2(au; + zFi) ) (6^) 



and where 

o2 



If ^ g ( ^^ f{E)-f{E + hu) f 1 l\ 



Here, we used 



A{E + rwj)A*{E) fiuj + iV \A{E + hu;) A*{E) 



(67) 



It turns out that g^{u!) is the ac-admittance for the symmetric barrier 
without gate (C = 0). This is is readily verified since for C = and for 
Oil = a22 it holds 

Giu.) ^ Gi, (C.) = ^e"?oV"e ^?re = ^H^) ■ (68) 

In the case of C = it is thus sufficient to discuss the symmetric barrier. 
We consider first this case and assume zero temperature. The integral over 
the difference of the Fermi functions reduces then to a simple integral from 
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^(e?) _ fii^ to n^^'^^ and can be carried out analytically. We find for the real 
part and the imaginary part of the admittance 



^ r M r / ,AE + nu^ ,AE-?iu;. 

Re[g = — —- I arctan( — — } — arctan( 



^ 4^0; V (r/2) ' ' (r/2) 



(69) 



T r // M r , /(|A|2-(na;)2)2 + (rna;)2 

where |A| = |A(/Lt(^*))|. Our results for g^{uj) are equivalent to earlier results 
of Fu and Dudley [6]. However, these authors neglect Coulomb interaction 
and, coincidentally, for a symmetric RTB, obtain the right expression only 
since they apply an antisymmetric bias zbAy/2 at the contacts. 

An expansion of Re[G(a;)] and Im[G(a;)] with respect to frequency yields 
the coefficients defined in Eq. (6) 

(71) 
(72) 

(27r|AI)2 ' 1^12 I • (73) 

We see, that E and K change sign. In particular, for the symmetric case 
this occurs if the value of the transmission probability T equals 1/2 and 
3/4. respectively. In Fig. 6a we plotted the real part and the imaginary part 
of the normalized admittance g\u))/g^ {Q) as a function of ^ti;/|A| for three 
different cases T = 1 (solid), T = 0.5625 (dashed), and T = 0.25 (dotted). 
For T < 3/4 there is a peak in the conductance which belongs to the 
excitation of the resonance by an energy quantum fiLO ~ |A[. We mention 
that the normalized admittance is determined by a single parameter 7 = 

r/iAi. 




We discuss now the effect of the capacitance C. In principle, by using 
the Eqs. (64) it is possible to derive expressions for the conductances (28). 
Here, we do not display the lengthy formulas but we rather present the 
results in a figure. A comparison of the admittance for (7 = (thin dotted 
curves) and for finite C (thick curves) is shown in Fig. 6b for the resonant 
case T = 1. One sees that for the small-capacitance case shown in the figure 
the Coulomb-coupling to the external gate enhances the capacitivc part of 
the admittance. Indeed, an expansion with respect to the capacitance yields 

GUu:) = i-ir+^Gico) - iuC^ + 0{C') . (74) 
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Figure 6. Frequency dependent admittance of the symmetric RTB. a) charge neutral 
(C = 0) case, 7 = r/|A|; b) Influence of a weak gate: Gn (solid) and G12 (dashed) for 
G|A|/e^ = 0.01, T = 1 (dotted: C = 0). 



In particular, one concludes J2af3 ^ai3 = —i^^C + O(C^) which follows also 
directly from the fact that the admittance '}2ai5 understood as 

a resistor (with an external admittance Y^^p ^a/s) series with a capacitor 
(with a capacitance C). 

The characteristic time-scales of tunneling problems [49] are a subject 
of considerable interest. With the analysis given above we are now able to 
identify characteristic times for the electric problem. We remark that the 
electrical problem investigated here leads to answers not for single electron 
motion, but for the collective charge dynamics. Consider the case of a RTB 
in the zero capacitance limit. From the expansion of the conductance up to 
quadratic order in frequency we can identify two characteristic frequencies 
or time-scales. There exists a frequency ooi at which the magnitude of the 
dc-current and the displacement current arc equal. The frequency LOi is 
thus determined by ri = \E\/G^, where E is the omittance and the 
dc-conductance. The time-scale ri is thus a generalization of the RC-time. 
For the RTB we find that this time is equal to ri = h/T at resoance 
and is also equal to ti = Ti/T far away from resonance. It vanishes for a 
Fermi energy at which the transmission of a symmetric barrier is equal 
to 1/2. A second characteristic frequency is obtained by comparing the 
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displacement current determined by the emittance E with the second order 
dissipative term K. The second characteristic time T2 is thus determined 
by T2 = At resonance this time is given by T2 = {l/3){h/T). This 

time tends to zero for Fermi energies far away from resonance. It is singular 
at the Fermi energies at which the emittance vanishes and it is zero at the 
Fermi energies for which K is zero. We conclude by mentioning that for the 
case in which the capacitance is not zero, the above consideration has to 
be applied to each element of the admittance matrix separately. In analogy 
to the tunneling time problem, where there exist characteristic times for 
traversal and reflection, from left and from right, the electrical problem will 
then be characterized by a set of characteristic frequencies or time scales 
for each admittance element, related only by sum rules due to current 
conservation and gauge invar iance. 

4. Conclusion 

In this work we presented a theory for the admittance and the nonlinear 
transport of open (i.e. connected to different reservoirs) mesoscopic con- 
ductors. We have emphasized the application of this theory to a number of 
systems of current interest, like quantum wires, quantum point contacts, 
and resonant tunneling barriers. Our emphasis has been to derive results, 
which even so they might not be realistic in detail, nevertheless capture the 
essential physics. Our results should be useful both for comparison with 
additional theoretical work and with experimental work. 

Due to the limitations of space, we have not reviewed a number of 
closely related subjects. The approach discussed here has been applied to 
ac-transport in two-dimensional electron gases in high magnetic fields [9, 38] 
which is particularly interesting in the regime where transport is dominated 
by carrier motion along edge states and by charging and de-charging of edge 
states. Of fundamental interest are Aharonov-Bohm effects in capacitance 
coefficients if one of the capacitor plates has the form of a ring [50] or for 
rings between capacitor plates [51]. As in dc-transport these interference 
effects are closely related to the sample specific fluctuations away from 
an average behavior. For a chaotic cavity capacitance fluctuations away 
from an ensemble averaged capacitance have been investigated by Gopar 
and Mello and one of the authors [52]. The ac-rcsponsc and its fluctua- 
tions of a chaotic cavity connected to two leads and coupled capacitively 
to a back gate has been investigated by Brouwer and one of the authors. 
This system is particularly interesting since the ensemble averaged quanti- 
ties exhibit weak localization corrections [53]. Since the weak-localization 
effect is not associated with a net charge accumulation the ac-response 
exhibits in addition to a Coulomb charge relaxation pole also a pole for un- 
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charged excitations. Scholler [54] has investigated the dynamic capacitance 
of a one-dimensional wire. The dynamic capacitance exhibits interesting 
structure due to the plasma-modes of the onc-dimcnsional wire. This list 
illustrates that there are many avenues to extend the work presented here. 
The application of the theory to more realistic models is a quite challenging 
undertaking but very likely also full of rewards. 
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